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Shaken optical lattices permit to coherently modify the tnnneling of particles in a controllable 
manner. We introduce a general relation between the geometry of shaken lattices and their ad¬ 
missible effective dynamics. Using three different examples, we illustrate the symmetries of the 
emerging tunneling rates. The results provide a clear framework to understand the relation between 
lattice geometry and accessible dynamics, and a tool to straightforwardly derive truncated effective 
Hamiltonians on arbitrary geometries. 


I. INTRODUCTION 

Advances in the control of quantum systems [1] have 
permitted to move the scope of experiments from corrob¬ 
orating fundamental effects to, most recently, designing 
systems that target specific effects such as geometric frus¬ 
tration [2, 3] or topological energy bands [4-7]. 

A common approach to modify the properties of a sys¬ 
tem is by means of external periodic driving [8, 9]. In a 
suitable fast-driving regime, driven systems can be used 
to mimic the dynamics of time-independent systems with 
properties determined by the amplitude and frequency of 
the driving [10-14] or, more generally, its temporal shape 
[15, 16]. 

In the domain of many-body physics, shaken optical 
lattices offer a particularly tunable platform for quan¬ 
tum simulations. Pioneer experiments in this context 
focused on the effects of the lowest order term in a high- 
frequency expansion of the effective Hamiltonan, and led 
to the observation of coherent suppression of tunneling 
[12] and synthetic magnetism [15, 17]. Most recently, 
the interplay of the leading and higher order terms has 
been exploited in order to simulate tunneling processes 
that lead to topological energy bands [5]. Taking advan¬ 
tage of such higher order terms broadens the possibilities 
for quantum simulations substantially, since the range of 
available physical processes that contribute to the effec¬ 
tive Hamiltonian grows with this order. 

The increasing possibilities, however, also come with 
increasing requirements towards the driving that needs 
to be chosen such that desired processes enter the effec¬ 
tive Hamiltonian with substantial weight whereas unde¬ 
sired processes are suppressed. This can be done with 
an appropriately chosen temporal shape of the driving 
[18], but the underlying geometry of a lattice provides 
limitations, and also potential for such control. 

In this article we investigate the influence of the lat¬ 
tice geometry on the existence or absence of higher or¬ 
der terms in the effective Hamiltonian. Specifically, we 
discuss how the symmetries of the lattice geometry lead 
to strict selection rules resulting from destructive inter¬ 
ference. These rules indicate which lattice geometries 
exclude the realization of certain processes and how un¬ 
desired processes can be avoided completely. In Secs. H 
and HI, we briefly review the notion of effective Hamilto¬ 


nians and the model of shaken lattices, respectively. The 
main results on the geometry dependence of the effective 
dynamics are presented in Sec. IV. Finally, in Sec. V 
we exemplify the symmetries of the emergent tunneling 
rates with three different geometries: a zig-zag chain, a 
kagome lattice and a Lieb lattice. 

II. EFFECTIVE HAMILTONIANS 

Effective Hamiltonians of driven systems naturally 
arise as consequence of Floquet theorem [19], which 
asserts that the time-evolution operator of a periodic 
Hamiltonian H(t) = H{t + T) can be expressed as 

U{t) = (0), (1) 

with a periodic unitary Upit) = Upit + T) and an Her- 
mitian time-independent operator i7eff- This operator 
defines the effective Hamiltonian of the system in the 
gauge Uf{Q) = In a suitable fast driving regime, 

where the driving frequency uj = 2 tt jT is the largest 
energy scale of the system, is nearly constant 

{Up{t) « Up(0)) with fast but small fluctuations. Any 
F satisfying limu,_>.oo F = 0 ensures that Up{t) approxi¬ 
mately equals that identity, such that e”*'^"*'*'* provides a 
good approximation of U(t). 

The representation in Eq. (1) can be found perturba- 
tively in powers of uj~^ using different methods [8, 20-22]. 
The effective Hamiltonian may then be expressed as 

+ ( 2 ) 

with the superscript indicating the order 0{uj~^) of 
each term. Similarly, the generator of the gauge can be 
written as 

F = -k -k , (3) 

where the expansion contains no zeroth order term F^^^ 
due to the requirement that Up remains close to the iden¬ 
tity. If a driven Hamiltonian leads to a zeroth order con¬ 
tribution F^^\ as in the case of strongly shaken lattices 
that we shall consider, one can work in a moving frame 
defined by a unitary Ui{t) where no such complication 
arises. 
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The two leading-order terms of the effective Hamilto¬ 
nian then read = Hq and 

oo 

= - E - ^o]) . ^ 

( 4 ) 

+ i[F^^\Ho], 

with the Fourier components H„ = 

(l/T)and the lowest-order term of 
the gauge generator The choice of gauge enters in 
the first-order expression of the effective Hamiltonian 
but it does not affect its eigenvalues up to the order of 
expansion. Typical conventions are the Floquet gauge 
F = 0, or the choice [8] 

fW ( 5 ) 

Lun 

n 

which leads to an effective Hamiltonian Fl^^g with no de¬ 
pendence on the static Fourier component. 


III. THE MODEL 


We consider non-interacting particles on a periodically 
driven lattice described by the Hamiltonian 

Hit) = H,+Hait). (6) 


The system Hamiltonian Fig characterizes the tunneling 
of the particles on a lattice with d sites in each unit cell 
(d-point basis). In a single-band tight-binding descrip¬ 
tion, it reads [10] 

F[g = ^ ^ c| Jjj Cj , (7) 

ij 


in terms of the dxd tunneling matrices Jij and the vector 
creation and annihilation operators 


= 

(c^ 

• ) 

1 

V 2,Si 5 

/ 

\T 

= 

(G,si ) ■ ■ 

) G,Sd) 


( 8 ) 


The indices i and j denote sites of the underlying Bra- 
vais lattice and the pair of indices {i, s^} characterize 
the physical lattice site Sk of the ith unit cell. The ma¬ 
trix elements Jij\spSg thus indicate the tunneling am¬ 
plitudes among physical sites. Infinite lattice size or 
periodic boundary conditions are assumed. The lat¬ 
tice site creation and annihilation operators satisfy 
the usual (anti-)commutation relations [ci^Sj,,Cjs ]± = 
Sij,s s ) where the commutator applies to bosons and the 
anticommutator to fermions. 

The time-dependent part of the Hamiltonian F[^{t) de¬ 
scribes the periodic driving of the system and is given by 
[23] 


HS)=Y.^\Emc^ ( 9 ) 

i 


with the diagonal matrix 

E,{t) = diag(-r,,si • r(t), • • • , • F{t)) (10) 

expressed in terms of the external driving force F(t) and 
the lattice site positions . We focus on a strong driv¬ 
ing regime where the amplitude of the driving force is 
comparable to the driving frequency. In order to con¬ 
fine the magnitude of the time-dependent oscillations it 
is convenient to perform a unitary transformation [24] 
and work in the frame defined by 

Ui{t) = exp ^c|0,(<)c,^ , (11) 

where Qi{t) = f* E^(T)dT - (f^ Ei(T)dT}T and (-jr de¬ 
notes time average. The Hamiltonian in Eq. (6) then 
transforms to E(t) = UiHUj — iUjdtUj, which we shall 
regard as our starting-point Hamiltonian for the remain¬ 
der of the article. It reads 

H{t) = Y^c\G,j{t)Cj ( 12 ) 

ij 

in terms of time-dependent matrices Gij{t) = 
g-jei(t) j_g 2 ej(t)_ Importantly, the Hamiltonian in Eq. 
(12) has the same translational invariance as the under¬ 
lying Bravais lattice of Fig- 


IV. GEOMETRY DEPENDENCE ON THE 
EFFECTIVE HAMITLONIAN 

The effective Hamiltonian of the shaken lattices de¬ 
pends crucially on the geometry of the lattice. In this 
section, we will explore this geometry dependence and 
put special emphasis on the symmetries of the geometry- 
dependent tunneling processes. 

With a high-frequency expansion of the effective 
Hamiltonian, as introduced in Eq. (2), the lowest or¬ 
der term corresponding to the Hamiltonian F[{t) in Eq. 
(12) reads 

(13) 

■ij 

where are the average or static Fourier 

components of the time-dependent matrices Gij{t) = 
Consequently, the tunneling rates of the 
leading order term in the effective Hamiltonian change 
with respect to the undriven ones and possibly describe 
different physics, but no new tunneling processes ap¬ 
pear. For instance, if the undriven system only contains 
nearest-neighbor (NN) tunneling, then the Hamiltonian 
necessarily contains only nearest-neighbor tunneling 
as well. That is, the leading order effective Hamiltonian 
cannot be used to engineer tunneling processes that are 
not present in the undriven system. 

New tunneling processes that are not present in Fig can 
nonetheless emerge in higher-order terms of the effective 
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(k) 

Hamiltonian , k > 1. We find that the existence of 
such higher order terms is closely related to the lattice 
geometry. 

(i) Bravais lattices. For Hamiltonians defined on a 
Bravais lattice with a one-point basis, the effective 
Hamiltonian is exactly given by H^ff = with 
all higher order terms necessarily vanishing = 
0, fc > 1, independently of the driving force. 

(ii) Non-Bravais lattices. For Hamiltonians defined 
on a Bravais lattice with a d-point basis, d > 1, 
higher order terms in the effective Hamiltonian ex¬ 
pansion can potentially appear depending on the 
specific time-dependence of the driving. 

This result, derived in Appendix A, shows a straight¬ 
forward distinction between those lattices geometries for 
which the effective Hamiltonian unavoidably has the 
same structure as the undriven one, and those for which 
the effective Hamiltonian can contain new processes and, 
thus, give rise to a more versatile platform for quantum 
simulations. 


Symmetries of the first-order term 

Given a non-Bravais lattice, the lowest order term of 
the effective Hamiltonian with emergent tunneling pro¬ 
cesses is given by H^g in Eq. (4). At this order, the par¬ 
ticular gauge defining ideff appears explicitly and needs 
to be specified. We choose it to be of the general form 

= (14) 

UJ 

where Fij are some tunneling matrices with the same 
translational symmetry as those of the undriven system, 
i.e. Jij. This includes, in particular, the Floquet gauge 
and the one defined by Eq. (5). For those gauges, the 
first-order term effective Hamiltonian can be written as 

(15) 

ij 

with the tunneling matrices Bij = Vuj and 

oo 

=E — 

n=l ^ (16) 



The emergent tunneling processes can be interpreted 
as the sum of all possible virtual two-step processes Vuj, 
where a particle tunnels from the jth unit cell to the zth 
unit cell through an intermediate unit cell 1. 

In order to inspect fundamental properties of these pro¬ 
cesses, we consider the matrix elements Bij\g^s^, which 
determine the tunneling amplitude for a particle at a 


physical lattice site jq = {j, s^} to tunnel to the site 
ip = {i, Sp}. They can be written as the sum of all pos¬ 
sible virtual two-step processes 

BijlspSg — 'y ( Pi-nlrjg J (1'^) 

where describes the tunneling amplitude of a par¬ 

ticle tunneling from the site jq to the ip via an inter¬ 
mediate site Ir = Explicit expressions for 

can be readily derived using Eq. (16) and are given in 
Appendix B. 

An important fundamental symmetry of the quantities 
Pipii- 3 q appears when the tunneling amplitudes of the un¬ 
driven system Hg only depend on the relative position 
between the sites, i.e. 

•/iilspS, = • (18) 

This condition is necessarily satisfied for Hamiltonians 
defined with a one-point basis (for which Fl^ vanishes) 
and for many other lattices such as the hexagonal. Even 
though some non-Bravais lattices like the kagome would 
admit a different configuration of the tunneling rates, the 
condition in Eq. (18) is both usually considered in the¬ 
oretical work [25-27] and also implemented experimen¬ 
tally with optical lattices [28]. When Eq. (18) applies, 
each of the virtual two-step processes can be character¬ 
ized by only two vectors {a^, a.j} independently of the ini¬ 
tial point. Consequently, the tunneling amplitudes 
in Eq. (17) can be written as a function 

Piplrjq =/3(3-i)Sj) (19) 

of the two vectors and a^ = Vi^g^ — 

The tunneling amplitudes I5{sii, SLj) then satisfy the 
symmetry 

/3(ai,aj) =-/3(aj,ai), (20) 

independently of the lattice geometry or external driving 
force, as demonstrated in Appendix B. Eq. (20) is di¬ 
rectly related to the fact that, for Hamiltonians defined 
with a one-point basis, Fljjg vanishes, as described in Fig. 
1. The symmetry in Eq. (20), however, also appears in 
Hamiltonians defined on non-Bravais lattices and leads 
to important constraints on the emergent tunneling pro¬ 
cesses. In particular, Eq. (20) leads to important selec¬ 
tion rules that depend on the specific lattice geometry, 
as exemplified in the next section. 

V. EXAMPLES 

Here we illustrate the symmetries introduced in Sec. 
IV for the gauge defined by Eq. (5) and for three dif¬ 
ferent non-Bravais geometries: a zig-zag chain, a kagome 
lattice and a Lieb lattice. We consider that the under¬ 
lying undriven system of all the examples contains only 
NN tunneling with isotropic rates satisfying the symme¬ 
try described in Eq. (18). 
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FIG. 1. Sketch of the triangular lattice and its two-step pro¬ 
cesses. All two-step processes can be defined either by only 
one vector {ai,ai} (red ©), or by a pair of vectors {ai,aj}, 
Ri yf Rj (colors green Q), purple (3) blue @) and its conjugate 
{aj,ai} (dashed colors). For the former case, Eq. (20) im¬ 
plies that the corresponding tunneling rate vanishes. For the 
latter, Eq. (20) implies that the two tunneling process of the 
pair interfere destructively. This yields a vanishing first-order 
effective Hamiltonian — 0, in agreement with the gen¬ 

eral geometrical distinction in the beginning of Sec. IV. This 
characterisitic is not specific to the triangular lattice but ap¬ 
plies to all lattices defined on a one-point basis. In particular, 
it is not restricted to nearest-neighbor tunneling. 


A. Zig-zag chain 


A zig-zag chain can be defined in terms of a one¬ 
dimensional Bravais lattice with primitive vector b, and 
two sites per unit cell denoted by Si and S 2 , as sketched 
in Fig. 2. For this geometry, the Hamiltonian H{t) in 
Eq. (12) becomes 


H(t) = 'Vc]; Go(t)c 


+ 




■^ri+b 


Gi{t)Cr^ + H.C.^ 


( 21 ) 


Go(t) = 


, Gi(t) = 


0 

0 0 


( 22 ) 


with the matrices 

0 g-RiitV 

v5ai(i) 0 

defined through the time-dependent tunneling rates 

(23) 

in terms of the tunneling rates jo, the time dependence 
Xk{t) = /q dr F(r) • aj, - ^ dt dr F(t) • a^, and the 
vectors a^ depicted Fig. 2. 

The lowest order term of the effective Hamiltonian 
Ffgg is readily given by the temporal average of Eq. (21) 
and reads 


4.^0Cr. + 


<+b^lCr 


+ H.C.) (24) 



FIG. 2. Sketch of a zig-zag lattice and its emerging pro¬ 
cesses. The vector b defines the primitive vector of the one¬ 
dimensional Bravais lattice. The vectors Ri and R 2 connect 
neighboring sites si and S 2 , which are indicated with red and 
green dots, respectively. The purple © arrows indicate the 
virtual processes giving rise to the on-site energy bo. Blue © 
arrows illustrate two emergent next-nearest-neighbor tunnel¬ 
ing with opposite rates. 


with the tunneling matrices 
^0 


Aq — 


0 

0 


— 


0 9l 

0 0 


(25) 


Consequently, contains renormalized NN tunneling 
rates that in general depend on the direction of tun¬ 
neling Rk- 

The first-order effective Hamiltonian term , calcu¬ 
lated from Eq. (15), becomes 

= E 4^0Cr. + E (4,+bSlCr. + H.C.) (26) 


with the matrices 
Bo = 


bo 0 
0 -bo 


Bi = 


bi 0 
0 -bi 


(27) 


The effective on-site energies and NNN tunneling rates 
are given by 


bo — + (3{—a2,^2) (28) 

6i =/3(ai,a2), (29) 

in terms of 

OO - 

/3(a„ a,) = ^ — ) . (30) 

TLUJ \ j j / 

n—1 


Eq. (30) characterizes the contribution to the effective 
rates arising from particles following the virtual process 
{a^, Rj}, i.e. tunneling from a site r to r + Ri + Rj through 
the intermediate site r -1- a^, as described in the previous 
section. 

The on-site energies and tunneling of the effective 
Hamiltonian have a clear interpretation. The on¬ 
site energy bo of sites Si, for instance, emerges as a con¬ 
sequence of the sum of the two processes {ai, — ai} and 
{—a 2 , a 2 }. On the other hand, the rate bi describing the 
tunneling between sites si of neighboring unit cells, re¬ 
sults from the process {ai,a 2 } only, as depicted in Fig. 
2. Similarly, the on-site energy —bo and rate —bi in 
Eq. (27) result from processes with the inverse order. 
The different signs thus arise from the general symme¬ 
try introduced in Eq. (20), in agreement with Eq. (30). 
The symmetries of the emergent tunneling imply that 
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Eq. (18) no longer holds for the tunneling of the effec¬ 
tive Hamiltonian, as the rate for particles to tunnel from 
a site r to r -I- ai -I- a 2 is different to the rate for particles 
to tunnel from r -|- ai to r -|- ai -|- a 2 . If ai = a 2 , a lin¬ 
ear Bravais lattice is recovered and Eq. (30) implies that 
bo = bi — 0 consistently with the geometrical distinction 
in the beginning of Section IV. 


B. Kagome lattice 


The kagome lattice is one of the most studied two- 
dimensional non-Bravais lattices, largely due to the cen¬ 
tral role it plays in the context of geometrically frustrated 
physics [25-27]. Recently, this geometry has been exper¬ 
imentally realized using ultra-cold atoms on an optical 
lattice [28], which suggests that its driving could be im¬ 
plemented in a near future. The kagome lattice consists 
of a triangular Bravais lattice with primitive vectors bi 
and b 2 , and a three-point basis, as depicted in Fig. 3. 
For this lattice, the Hamiltonian in Eq. (12) becomes 

H{t) =Y^cl^Go{t)Cr^ 

* 3 (31) 

+ J2Y1 (4.+b,G'fe(t)Cr, -b H.c. j 

i k—1 


with the vector ba = —bi — b 2 , the matrices 


/ 0 5a3 \ /O 0 0 \ 

Go{t) = 5ai 0 5_a2 , Gi(t) = 0 0 5a2 

\9-.s 9^. 0 J yo 0 0 / 

/ 0 0 0 \ /O 0 \ 

G2{t) = 0 0 0 , Goit) = 0 0 0 

Vffaa 0 0/ yO 0 Oj 


and the same time-dependent tunneling rates (t) intro¬ 
duced in Eq. (23). The vectors connecting neighboring 
sites are represented in Fig. 3. 

The leading effective Hamiltonian is readily ob¬ 
tained from the time average of Eq. (31) and does not 
contain new tunneling processes. 

The first-order term , on the other hand, is given 
by 


=E4^0Cr. 


-b 


EE(d 


Bi-Cj. + H.c. I . 




i k—1 


The tunneling matrices are 

/O bl bo 
Ho = hi 0 bl 
\bl b2 0 

0 0 0 

Uf* f 


\ J 

( 0 

0 - 

-b'o*\ 

, Bi=\ 


0 - 

-b2 



0 

0 J 

\ 1 


-bl 

b'A 

, Bo =\ 

0 

0 



^0 

— h'* 

0 / 


(33) 


(34) 



FIG. 3. Sketch of a kagome lattice, bi and b2 denote the 
two primitive vectors of the underlying triangular Bravais lat¬ 
tice. Each unit cell contains three sites Si, Si and sa in red, 
green and blue, respectively. The vectors a*, connect nearest- 
neighbor sites. The virtual paths of two emergent nearest- 
neighbour tunneling with opposite rates are indicated with 
green 0 arrows. 


defined in terms of the effective NN tunneling rates 


6 1 = ^(-a3,-a2), (35) 

62 =/3(-ai,-a3), (36) 

63 =/ 3 (-a 2 ,-ai), (37) 

the NNN tunneling rates 

61 =/3(-a3,a2), (38) 

62 =/3(-ai,a3), (39) 

^3 =/3{-a2,ai), (40) 


and with /3{a.i, aj) in Eq. (30). As depicted in Fig. 3, the 
relative signs in Eq. (34) are a particular manifestation 
of the fundamental symmetry in Eq. (20) and appear 
due to the virtual path that particles follow. 

There are two basic differences between the emergent 
processes of the zig-zag chain and those of the kagome 
lattice. On the one hand, new NN tunneling appears in 
the kagome lattice as a consequence of two consecutive 
virtual NN tunneling in a triangular geometry. On the 
other, no on-site energies appear because the different 
contributions interfere destructively. For each process 
{a^,—a^} contributing to the on-site energy, there is a 
conjugate process {—ai,ai}. This exemplifies well how 
different non-Bravais lattices can induce different effec¬ 
tive tunneling processes due to the specific details of the 
geometry, which opens the possibility to designing spe¬ 
cific geometries that yield desired tunneling processes by 
driving them. 
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that can lead to a finite contribution of the emerging tunnel¬ 
ing in the effective Hamiltonian correspond to all pos¬ 
sible combinations of ±ai and ±a 2 . Virtual paths of two 
next-nearest-neighbor tunneling are indicated with blue 0 
and orange (2) arrows. 

C. Lieb lattice 

With the previous examples, we have demonstrated 
a clear interpretation of the tunneling rates of the first- 
order effective Hamiltonian in terms of particles virtually 
tunneling through an intermediate site. In fact, this ge¬ 
ometrical interpretation permits to infer the tunneling 
rates of effective Hamiltonians of shaken lattices essen¬ 
tially by only inspecting the geometry of the lattice, as we 
will describe here with particles on a shaken Lieb lattice. 

The Lieb lattice is a two-dimensional face-centered 
square lattice, defined in terms of a square Bravais lattice 
and a three-point basis, as depicted in Fig. 4. Analo¬ 
gously to the previous sections, the driven Hamiltonian 
H{t) in Eq. (12) is defined by its time-dependent tunnel¬ 
ing rates ga.i(t) in Eq- (23), characterizing the tunneling 
of particles from r to r -b a^. 

According to Eqs. (17), (19), and (30), the effective 
tunneling rate for particles to tunnel between two ar¬ 
bitrary sites can be determined one by one by inspect¬ 
ing all the possible two-step processes connecting them 
and adding up all the contributions. That is the rate 
6(r,r -b Hi -b aj) for particles to tunnel (or stay) from a 
site at r to r -b ai -b aj is given by 

b(r,r + a, + aj)= ^ ^(a^aj) (41) 

{ai.aj} 

where the sum is performed over all possible two-step 
processes and, in the gauge defined by Eq. (5), /3(ai, aj) 


is given by Eq. (30). The symmetry in Eq. (20) essen¬ 
tially establishes which processes necessarily vanish and 
which processes can be finite for specific driving forces. 
For the Lieb lattice it can be thus inferred that only 
NNN tunneling between green and blue sites can ap¬ 
pear in the effective Hamiltonian as described in 

Fig. 4. Moreover, the corresponding tunneling rates are 
,5(ai,a2), /3(a2,ai), /3(-ai,a2), ,5(a2,-ai) and its com¬ 
plex conjugate values. All other processes such as on¬ 
site energies or NNN tunneling between translationally 
equivalent sites (i.e. between sites with the same color in 
Fig. 4) do not appear because the contributions either 
are zero or sum up to zero, as analogously described in 
Fig. 1. 


VI. CONCLUSIONS AND OUTLOOK 

The structural distinction of effective Hamiltonians of 
shaken lattices in terms of the lattice geometry described 
in Sec. IV, establishes a general framework to understand 
the central role that the geometry of shaken lattices plays 
on its effective dynamics. 

The geometry dependence of the effective tunneling 
rates can be interpreted in terms of path-dependent vir¬ 
tual tunneling processes that interfere with each other. 
Depending on the details of the geometry, this interfer¬ 
ence can then lead to vanishing or finite effective tunnel¬ 
ing rates. The resulting selection rules open new possi¬ 
bilities for quantum simulations by providing a tool to 
readily identify lattice geometries that contain only de¬ 
sired tunneling processes. 

In a fast-driving regime, the lowest-order contribution 
to the effective Hamiltonian with geometry-dependent 
tunneling is given by the first-order term of a high- 
frequency expansion. At this order, the emerging tun¬ 
neling processes result from two consecutive virtual tun¬ 
neling processes of the undriven system. Consequently, if 
the undriven system contains only nearest-neighbor tun¬ 
neling, the possible emerging processes are on-site ener¬ 
gies and next-nearest-neighbor tunneling. The geometry 
dependence of the aforementioned next-nearest-neighbor 
tunneling can then be employed to engineer topological 
energy bands with non-zero Chern number [5, 18, 29, 30]. 
Our general characterization thus permits to identify a 
wide range of lattice geometries that can host topological 
phases. 

Highly interesting for actual quantum simulations is 
the inclusion of interactions or spin-dependent tunnelling 
processes [5, 31, 32]. The identification of structure- 
dynamics relations established here is rather general and 
not limited to the specific case of non-interacting spin¬ 
less particles. Since the approach is based on general 
geometric and algebraic properties extensions to systems 
with different underlying dynamics may readily be ob¬ 
tained and also extensions to higher order processes are 
rather straightforward. The possibility to identify struc¬ 
tural properties of a quantum simulator based on its un- 
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derlying geometry is thus expected to aid substantially 
in the design of optimal implementations thereof. 
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Appendix A: Geometry distinction 

The geometrical distinction established in Section IV is 
based on the commutativity of the Hamiltonian at differ¬ 
ent times and can be derived with the Magnus expansion 
[20, 33]. 

The Magnus expansion approximates the time- 
evolution operator as U(t) = with 

M{t) = Moit) + Mi(t) -!-•••, (Al) 

where MQ{t) = J* dtiH{ti) and higher-order contribu¬ 
tions Mk(t) are given in terms of fc — 1 time integrals 
of /c-fold commutators. The effective Hamiltonian of 
the driven system can be then obtained order by order 
through 

HeB = ^uI{0)M{T)Uf{0) (A2) 

for any chosen gauge Uf{0)- 

A sufficient condition for all Magnus terms Mk{t), with 
fc > 1, to vanish is that [H{ti), H{t 2 )] = 0 [33]. Conse¬ 
quently, 

[H{h),Hit 2 )] = 0 => ileff = (A3) 

with all higher order terms of the effective Hamiltonian 
exactly vanishing. 


Using the (anti)commutation relations of the creation 
and annihilation operators of the Hamiltonian in Eq. 
(12), the commutator of the Hamiltonian can be writ¬ 
ten as 

[H{ti),H{t2)] = J2clMijiti,t2)cj, (A4) 

with 

M,jih,t2) = Y, iGuih)Gij{t2) - Guit2)Gi,{h)lA5) 

i 

The translational symmetry of the Hamiltonian permits 
us to use the identity 

'Y^GiiDij — DaCij) = 'Y[Gii,Dij\, (A6) 

i i 

which holds for translationally invariant matrices and 

Dij as derived in Appendix C. With this relation, Eq. 
(A5) van be rewritten in the more convenient form 

= Y\GaiG),Gi,{t2)l (A7) 

... * 

from which the distinction between Bravais and non- 
Bravais lattice follows: if the Hamiltonian is defined on a 
one-point basis, the Gij{t) are scalars and the commuta¬ 
tor necessarily vanishes. If it is defined in a d-point basis 
with d > 1, the Gij(t) are matrices and the commutator 
can be non-vanishing. 

This distinction considers the lattice geometry through 
the commutativity of the tunneling matrices but it does 
not take into account the specific time dependence of the 
Hamiltonian. It is thus possible that systems within the 
category (ii) in Sec. IV still satisfy ideff = due to 
particular driving profile. 

Appendix B: Symmetry of the tunneling rates 

Expressions for Bi i can be readily derived using 
Eqs. (16), (17) and (A6)! They read 




-g; 


n—1 

oo 


|spSrG;jjs,.Sg 


+ E - Gu)KsGlUs, - G° |.^.^(G-" - Gf^)UsJ 


(HI) 


+ ^{Fu\spSr.Glj\s^Sg Ggl\spSr.Flj\sr.Sg)- 


r 


plies for both Fij and G"-, i.e. 


FijlspSg — /ri,sp-rj,a, 

nn I _ n 

^iilSpSq — y : 


ri.s^-rT 


The symmetry P{a 4 ,aj) = —/3(aj,ai) introduced in 
Eq. (20) can be derived from Eqs. (Bl) and (18). The 
condition in Eq. (18) implies that the same relation ap¬ 


(B2) 

(B3) 




Eq. (B2) is a consequence of our choice of gauge, which 
has the same symmetries as J^. Eq. (B3), on the other 
hand, can be derived by calculating the matrix elements 
of Gij{t) in Eq. (12). 

Consider now two tranlstionally invariant matrices Cij 
and Dij with analogous symmetries as in Eqs. (B2) and 
(B3). We then define the quantity K(ai,a 2 ) as 

A^(ai,a2) Gil\s^s^]Dlj\s^s^ ^i[\spSr.Glj\s^Sq 

~ — driap-ria,,Cr,,^-rj(B5) 

— Ca2 dai da2 Cai 5 (B6) 

which explicitly satisfies K(ai,a 2 ) = —K(a 2 ,ai). Com¬ 
bining this result with Eqs. (Bl), (B2) and (B3), the 
symmetry /3{aLi,aj) = —/3(aj, a^) introduced in Eq. (20) 
directly follows. 


Appendix C: Identity for translationally invariant 
matrices 

Eq. (Cl) can be proven by exploiting the transla¬ 
tionally invariance of the Bravais lattice and reordering 
the terms of the sum. For convenience, we introduce a 
slightly different notation and denote translatinally in¬ 
variant matrices Cij by C(ri — r;), where the argument 
Ti — ri does not indicate an explicit dependence on the 
vectors but rather a label. With this notation, Eq. (A6) 


can be expressed as 

^(^(r, - ri)D{ri - r^) - D{ri - r/)C(r, - r^-)) = 

I 

J2[Ciri-ri),D{ri-rj)]. (Cl) 

i 

The translational symmetry of the Bravais lattice im¬ 
plies that, if the vectors and r_, denote the position of 
two Bravais sites, then for each Bravais site I at a posi¬ 
tion r; = r.i -I-Rh, there exists a “conjugate” Bravais site 
I* at a position r/. = — R^. Consequently R/^ = R^-;*. 

Moreover r;. = r; if and only if R^ = (r^ — ri)/2. 
Consider a term 

Ti = C{r^ - ri)D{ri - rj) - D{vi - Yi)C{vi - r^](C2) 

in the sum of Eq. (Cl) with fixed i and j. Then, either 
R/i = {vj — ri)/2 or R/^ ^ (r^ — ri)/2. The first case 
implies that the conjugate site of I is itself and thus T; 
can directly be written as the commutator 

Ti = [Cir,-ri),Diri-r,)], (C3) 

since R/i = Rji. The second case implies that r; ^ r;* 
and thus there there exists another term Ti* in the sum 
such that 

Ti + Tf = 

[C{vi - vi),D{vi - Tj)] -b [C{vi - vi*),D{vi* - rj)(]C4) 

Therefore, the entire sum can be rewritten as 

Y,Ti = Y,[C{v,-vi),D{vi-v,% (C5) 

i i 

which completes the proof of Eq. (Cl) and hence (A6). 
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